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Abstract 

We consider the scattering transform for the first order system in the plane, 




We show that the scattering map is Lipschitz continuous on a neighborhood of 
zero in L 2 . 

This paper gives an estimate for the scattering map associated to a first-order 
system 

Dip — Qijj = (1) 
in the plane. Here, D and Q are denned by 

D = ( 9 n a I and Q = 




and dx and d x are the standard derivatives with respect to x = x 1 + ix 2 and x. 
The entries of the matrix Q, g 1 (a;) and q 2 (x) are complex valued functions on the 
complex plane. (We will consistently use superscripts to indicate components.) The 
system ((TJ) was studied by Beals and Coifman [21 E] , and a number of other authors 
(see Fokas and Ablowitz [7j for an earlier formal treatment, Sung El ESj for a 
detailed rigorous treatment and the review articles j3J E] for additional references). 

* Author supported in part by the National Science Foundation, Division of Mathematical Sci- 
ences. 
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One beautiful feature of the scattering theory as developed by Beals and Coifman is 
the existence of a non-linear Plancherel identity relating the scattering data S, which 
is defined below, to the potential Q in ((IJ), 



S\ 2 d/i= / \Qrdfi. (2) 
c Jc 

Here, we use dfi to denote Lebesgue measure and 2 x 2-matrices are normed by 
1^1 2 = E 2 , fc= i \F jk \ 2 . The identity is valid for nice potentials Q satisfying one of 
the symmetries Q* = ±Q. Since the map Q — > S is non-linear, the identity (j2J) does 
not immediately imply the continuity of this map. The goal of this paper is to provide 
a proof of the continuity in the L 2 -metric of the map Q — > S for small potentials. 
Our methods apply equally to the inverse map S —> Q. We also compute explicitly 
the size of potentials for which our method works and obtain results for potentials 
with llQHa < y/2. 

Before describing the scattering theory and the proofs of our estimates, we indicate 
several reasons for interest in this system. The scattering theory for this system can 
be used to transform the Davey-Stewartson II system to a linear evolution equation. 
Thus, our estimates can be used to assert that for small initial data, the solution u 
depends continuously on the initial data in the L 2 -norm. However, it is not clear that 
this method produces a solution for all initial data in a neighborhood of zero in L 2 . 
The application to the Davey-Stewartson system seems to have motivated the work 
of Beals and Coifman, Fokas and Ablowitz and Sung. 

This system also appears in the work of the author and G. Uhlmann [3] who study 
the inverse conductivity problem in two dimensions and prove that the coefficient is 
uniquely determined by the Dirichlet to Neumann map for conductivities which have 
one derivative in L p , p > 2. This is the least restrictive regularity assumption on 
the coefficient which is known to imply a uniqueness theorem for the coefficient. The 
argument depends in a crucial way on the non-linear Plancherel identity (J2J). More 
recently, Barcelo, Barcelo and Ruiz |L have shown that the coefficient depends con- 
tinuously on the Dirichlet to Neumann map under the a priori regularity hypothesis 
that the coefficient has a gradient which is Holder continuous. One step of their 
argument is the observation that the scattering map is continuous if the potential is 
assumed to be in C e and compactly supported. The research reported here is perhaps 
a step towards relaxing this assumption. 

We begin the formal development by reviewing the notation of Beals and Coif- 
man's paper jH] . We let x = x 1 + ix 2 and z = z 1 + iz 2 denote variables in the complex 
plane. A family of solutions of the free system, Dipo = 0, is given by 

ipo(x,z) = 

We look for solutions of (JTJ) of the form ip(x, z) = rn(x, z)ip (x, z) where m approaches 
the 2x2 identity matrix as \x\ — > oo. A computation shows that m must satisfy the 
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equation 

D z m = Qm (3) 

where D z is the operator E~ l DE z and for z e C, E z is the map on 2 x 2 matrix- valued 
functions defined by 

E z f = r + A- Z i r. 

Here and below, we use / d and f° s to denote the diagonal and off-diagonal parts of 
the matrix / and A(x, z) = A z (x) is the matrix given by 



A(x, z) 



exp(ixz + ixz) \ _ / a 1 (x, z) 

exp(— ixz — ixz) J \ a 2 (x, z) 



The solution to (JHJ) is found by solving the integral equation 

m = l + G z (Qm). (4) 

In this equation and below, we use 1 to denote the 2x2 identity matrix. We use 
G z = E~ l GE z to denote the inverse of D z , where 

In our situation, with Q small, we can write the solution as a Neumann series, m = 
Y^k=o(G z Q) k (X) ■ Here and throughout this paper, we use Q to denote both the 
function Q and multiplication operator / — ► Qf. 

Continuing with the scattering theory, we can differentiate the solution with re- 
spect to z and obtain the <9 2 -equation 

d z m(x, z) = (Tm)(x, z) (5) 

where the right-hand side Tm is defined by 

Tm(x, z) = m(x, z)S(z)A(x, —z). (6) 

In the definition of T, S is the scattering data which is defined by 

S{z) = --J ( E z (Q(x)m{x, z)) d/j,(x) (7) 



7T 

and Jf = 2(J/) ofi with J the 2x2 matrix, 

J 2 V i ) ■ 

The scattering data S defined in (JJJ) is the function which appears in the left-hand 
side of (J2J). The identity (J2J) is proven for nice functions by Beals and Coifman in 0. 
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A detailed argument is given by Sung in El E>. for potentials Q which are in 
L 1 fl L°° and much of his work extends to potentials in L p H L p (p ^ 2). 

Our main result is to show that for Q small in L 2 , (and say Q G 5(C), the 
Schwartz space), the map Q — > S^Q) is Lipschitz continuous in L 2 and that a similar 
result holds for the inverse map S — > (5(5*). These results allow us to extend the 
maps to a neighborhood of zero in L 2 . In fact, the scattering map Q — > S(Q) and the 
inverse map S — > Q(S) are of the same form. (The more adroit normalization used 
by Sung makes this clear.) Thus, we concentrate most of our energy on the scattering 
map. 

We remark that the definition of S given in does not make sense if Q G L 2 , 
since I do not know if we can construct m without the assumption of some extra 
decay on Q. One assumption that will suffice is that Q G L p for some p < 2. 

Our proof is straightforward. We write m, the solution of (@J) as a Neumann series 
and substitute this series into the definition of S, (J7J). This gives a series of multilinear 
expressions in Q. We estimate each term and show that we can sum the series. The 
proof of these estimates depends only on the Hardy-Littlewood-Sobolev theorem, but 
requires a certain amount of persistence. 

The main ingredient is the following well-known result on the first order Riesz 
potential defined by 

Jc \x — y\ 

We will only consider R acting on scalar functions. 

Lemma 1 Let r and p satisfy r > and 1 < p/r < 2, then the map f — > [i?(|/| r )] 1 / r ' 
satisfies 

\\[R(\f\ r )] 1/r h<^r/p)^\\f\\ p 

where 1/p — 1/p = l/(2r). The constant a{9) is bounded for 9 in compact subsets of 
the interval (1/2, 1) and we have 

a (3/4) = 7i and hmsupa(#) < tv. 

6^3/4 

Proof. According to the Hardy-Littlewood-Sobolev theorem, on fractional integra- 
tion, (see [12] . for example), we have R : L s — ► L s where 1/s — 1/s = 1/2, provided 
1 < s < 2. We let a (1/s) denote the operator norm of this map on scalar- valued 
functions. 

Now, given p and r as in the Lemma, we define p by r/p — r/p = 1/2. We apply 
the Hardy-Littlewood-Sobolev theorem with exponents p/r and p/r to conclude that 

^J c [R(\fn(x)f r d^x)J /P < a{r/p) Qf \f{x)\ p d^x^ . 

Taking the rth root gives the first inequality in the Lemma. 

The operator norm for 1/s = 3/4 was computed by Lieb [S]. The behavior near 
8 = 3/4 asserted in the Lemma follows from Lieb's result, the boundedness for 1 < 
p < 2 and the Riesz-Thorin interpolation theorem. I 
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We now state the main result. 
Theorem 2 Let 

N = L 2 (C) n {F : \\F\\ 2 < V2 and F d = 0} 

then the maps S — > Q and Q — > S, defined initially on S(C)nN, extend continuously 
to N. In addition, 

S o Q = Id and Q o S = Id 
provided Q and S(Q) lie in N (respectively, S and Q(S) lie in N). If Q = ±Q* , then 



[ \Q\ 2 dfi = [ \S\ 2 d^. 
Jc Jc 



The proof will depend on multi-linear estimates involving repeated fractional in- 
tegration. For the study of these multi-linear expressions, we will work in a sequence 
of L p spaces using the exponents, Pj, Sj and r, defined by: 



p = 2 and s = 4. 



Then we define 



1 


4 







%' 


J 


1 


1 


1 


Pj+l 


Pj 


27? 


1 


1 


1 










s j 


2r~/ 


1 


1 


1 


§ 3 


s j 


2 



0,1,2,... 



These sequences satisfy 



Vj _ = 4 

Tj 3' 

1 1 3 

Pj s j 4' 

£i 4 

r) ^ 3' 



j = 0,l,2,... (9) 

j = 0,l,2... (10) 

j = l,2,... (11) 

j = 0,1,2,... (12) 

j = 0,l,2,... (13) 



.7 = 0,1,2,... (14) 
8j < 2, J = 1,2,... (15) 



Here and below, we use r' = r/(r — 1) to denote the dual exponent. The statement 
fll2|) is immediate from (jHJ and (j!3)l follows from the definitions of pj and Sj. The 
third, (JHJ), follows since (fTB*)) implies Sj = 4pj/(3pj — 4) and 1/H = (3pj — 4)/(3p 3 -). 
The last observation ()15|) follows since si = 12/7 and decreases. 
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The multi-linear expression we will study is Ik(t, qo, ■ ■ ■ q2k) which is defined for 
fc = l,2,...by 



h(t, qo, ■ ■ ■ , q2k) 

t(x - xi + x 2 - . . . - x 2k -i + x 2k ) IIj=o gjfoj) 
c 2fc +! b - xi\\xx - x 2 \ ■ ■ ■ b 2 fc-i - x 2k \ 



d/j,(x 0} . . .x 2k ). 



Also, we will use Io(t,q ) = J c t(x)qo(x) dfj,(x). The functions t and qj are scalar- 
valued. The main estimate for this expression is in the following Lemma. 

Lemma 3 For every e > 0, there exists a constant C e so that 

2k 

h{t, qo,---, q 2k ) < C e 7r 2k (l + e) 2k \\t\\ 2 ]J \\ qj \\ 2 . 

Before, beginning the proof of Lemma El we state and prove the following result. 
The proof of Lemma H3 amounts to applying this result k times. 

Lemma 4 Let j > 0, k > 1 and suppose that pj, Sj, s~j + i and r 3 - are exponents as 
defined in For t and qj non-negative functions on C, we have 

I k (t, q ,.-., q2k) < h-i{h, <?2<?2, 93, • • • , to) 
where the functions t\ and q 2 satisfy 

\\ti\\ Pj+1 < a(3/4)^i|| f || w 

||?2||« J+ i < «(3/4) 1/ ^a(s^ 1 )||gi|| 2 ||goL J - 

Proof. We consider the integral with respect to Xq in I k and apply the Holder in- 
equality using exponents rj and to obtain 

t(x - x\ + x 2 - . . . + x 2k )q Q (x ) . 



C \Xq — Xi 

< 



C \Xq — X\\ J 



l/r'. 
1 j 



x ( [ trS{x °- X ] + X2 ---- + X2k) d»(x )) 
\Jc bo — Xi\ J 

qi(xi) h(x 2 -x 3 + x 4 - ... + x 2k ). 



1/r, 



Here, the function t\ is defined by 

d/j,(w) \ 

\w\ 



h(x) = [R(fi)} l/r i (x) =( f f ^ X . , W ^ d/j(w 



One may make the change of variables w = x\ — xq in the integral with respect to xq to 

obtain this representation of ti. The function qi is defined by q\{xi) = [^(g^)] 1 ^ (xi). 
Now, we can rewrite the integral with respect to Xi in I k as 



q 1 [xi)q 1 [x 1/ 

\X2 — X\\ 



dn{xi) = R(q 1 qi)(x 2 ) = R{q 1 [R{q Q 3 )] 1/r 'i){x 2 ) = q 2 {x 2 ) 



which defines q 2 . Inserting the definitions of ti and q 2 into I k gives the inequality 
relating I k and Ik-i- 

We now establish the estimates for t\ and q 2 . The estimate 

¥i\\ Pj+i <<?>/ ±) x,r m\ Pj 

follows from LemmaQwith p = Pj and r = rj. The hypothesis rj/pj G (1/2, 1) follows 
from (|12j). To estimate q 2 , observe that 4/3 < Sj + i < 2 from (|T3|) and (|15p and thus 
we may apply Lemma ^ with p = Sj+\ and r = 1 to obtain the first inequality below. 
Next, we use Holder's inequality and Lemma ^ with p = Sj and r = r'- to obtain the 
second and third inequalities: 



\\R(qi[R(qh 1/r '0k +1 < a(«7i)ll«i[^?)] 1/rJ l 



< a^llftllall^)] 1 ^!!^! 

< «( S 7^ 1 )c(3/4) 1 /^ ||^! lUH^olU,- 

These are the estimates of the Lemma. I 

We are now ready to give the proof of Lemma El 

Proof of Lemma^ It suffices to prove the Lemma when all of the functions are 
positive. We first claim that if s = 4 or 4/3, we have 

2fc-l 

\I k (t,q ,...,q 2k )\ ^ Ce7r 2fc (l + e) 2fe ||t|| 2 || 9o || s ||?2fe|| s ' II Nk (16) 

3=1 

We use the sequence of exponents defined in (IH HTTf . By Lemma |U applied k times, 
we obtain 

hit, q , . . . , q 2k ) < Io(tk, q-2.kQ.2k) 
where the function t k satisfies 



\tk\\ Pk <a(3/^=^)\\t\\ Po 



and the function q 2k satisfies 



felk < a(3/4) E ^ [Haft 1 )] II WlMllo 



k-i j_ I k \ 2k-\ 

SO' 



\3=1 / 3=1 
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We write out the integral defining Iq and obtain 



c 

We have tk G L Pft , g 2 fc G L s ° and g2fc G L Sk where the exponents pk, s~k and s' satisfy 

111 1111 

— + — + - = — + -- - + - = 1. 

Pk s k s p s 2 s 

Thus Holder's inequality and the estimates for tk and q 2 k imply that 

k 2fc-l 

\h(t,q , ■ ■ ■ ,g 2 fc)| < a(3/4) fc J| ^(sJ^llgolUllfelU/sllilh II ll^'lh- 

i=i i=i 

Since — > 4/3 as j — > oo, Lemma [T] tells us that limsup^^ ai^sj 1 ) < ir. Thus, if 
e > 0, then we obtain (fTT)|) with s = 4. To obtain the result with s = 4/3, we simply 
need to change variables to replace by x 2 h-j- 

Now that we have the estimate we apply a simple interpolation argument to 
obtain the Lemma. We define an operator T by 

rr ( x /" i(^o -zi + --- + z2jfc)gi(ai)---g2&(a2A!) , / s 

Tq 2k (x ) = ; ; ; ; d/x(xi,...,x 2 fc). 

Jc2fe -^l| • • • \X2k-l -X 2 k\ 

According to the estimate pfij) . T maps L 4 — > L 4 and L 4//3 — > L 4 / 3 . Hence, by the 
Riesz-Thorin interpolation theorem, T is bounded on L 2 , with the same bound. This 
implies the Lemma. I 

Proof of Theorem 2. We begin by indicating how we construct the solutions m when 
the potential Q is nice, Q G 5(C), say, and \\Q\\ 2 < V2. We can write m as the infinite 
sum 

oo 

m = l + E(G.g) i (l)- (17) 
i=i 

We consider one entry in the matrix (G z Qy(l) with j = 2k. Note that for j even, 
the off-diagonal part (G z Q)i(l) is zero. Also, note that the operator E z involves 
multiplying by exponentials of modulus 1, and thus the entries of E z f have the same 
L p -norm as the entries of /. We apply Lemma Q with p = 4/3 and r = 1 to obtain 
that 

\\((G z Q) 2k )(l) n \U < ir- 2k \\R(Q 12 (...R(Q 21 )...))\U 

< (||Q 1 l2||Q 2 1 2 )*- 1 ||Q 12 ||2||Q 21 ||4/ 3 



< 



1 



illQllf^HQll 



4/3- 



The last inequality uses the elementary inequality that ||<5 12 ||2||Q 21 ||2 < |(IIQ 12 ||i + 
\\Q 12 1| 1). A similar argument holds for the remaining entries. Thus the infinite sum 



8 



in ()17|) will converge in L 4 if HQH2 < \/2. We substitute the sum for m into (J7J). 
Since Q and S are off-diagonal, only the even terms of the series for m are needed in 
the expression for S. Also, note that the operator G z = G when acting on diagonal 
matrices. Using these observations and that Q is off-diagonal, we can write 

— 2 r 00 
S(z) = — J A(x,-z)Q(x)y2(GQG z QY(l)dn(x) 

7T JC j=0 

r\ OO 

The term with j = in this series is essentially the Fourier transform. If we 
consider one entry, we obtain 



S$*(z)= / exp(-2i(a:V + x 2 z 2 ))Q 12 {x)dfi{x) 
and we have a similar expression for S 21 . Thus, 

— / \S \ 2 dfi = [ \Q\ 2 dfi 

TX Z JC JC 

by the standard Plancherel identity. The higher order terms require a bit more work. 
In order to simplify the notation, we consider the upper-right entry in S k 

Sk 2 ( z ) = ~9k f a 1 (-XQ+x 1 -x 2 +x 3 ... + X2k-i-x 2 k,z) 
n 2h Jc 2fc + 1 

Q 12 (x )Q 2 \ Xl ) . . .Q 21 (x 2k -i)Q 12 (x 2k ) 
x — — - — = dfi{x , . . . , x 2k ). 

(Xq - Xi)(Xi - X 2 ) ... (x 2 fe-2 - X 2 fc-l)(x 2 fc-l - X 2k ) 

We take a sufficiently regular (scalar-valued) test function T and write out the ex- 
pression for Si 2 to obtain 

J^T(z)Sl 2 (z)d{i(z) = J c2k+2 T(z)a 1 (-x + x\ - x 2 + x 3 . . . + x 2 k-i ~ x 2k ,z) 

Q 12 (x )Q 21 (x 1 )...Q 21 (x 2k ^ 1 )Q 12 (x 2k ) 
x — — r — — afi{z, x , ... , x 2k ). 

{Xq ~ -X 2 ) . . . (X 2fc _ 2 - X 2 fc-l)(^2fc-l - %2k) 

Note that if each of the functions T and Q are in L 1 D L°°, say, then the integral 
on the right of the previous expression converges absolutely. Thus, we may use the 
Fubini theorem to carry out the integration with respect to z first and obtain the 
expression 



1 



T(2(x -X1+X2- x 2k -i + x 2k )) 



Q 12 (x )Q 21 (x 1 )...Q 21 (x 2k . 1 )Q 12 (x 2k ) , . . 

Xtt —^7 r TZ Z w r dfiyXf), . . . ,x 2k ). 

{Xq ~ Xi){Xi -X 2 ) ... [X 2k _ 2 - X 2fc _i)(x 2fc _i - X 2k ) 
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Here, T is the Fourier transform of T and is defined by T(z) = J f(x)e dpi(x 
J f(x)a l (—\x, z) d/i(z). Thus we have that 



T(z)S 12 (z) dfi(z) 



<— I k (\T(2.)\,\Q 1 \\Q*\,...,\Q 12 \). 
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Now, Lemma El implies for each e > 0, there is a constant C e so that 

n -+- f) 2k 
l \s k h.<a K Z, } \\Q\\l k 



and hence the series 



2 A: 



oo 
fc=0 



converges in L 2 if \\QW2 < V%- 

We now turn to the proof that the map Q —>■ S is Lipschitz continuous. Let Q 
and Q be two matrix potentials. We write S = S(Q), then we have 

T(z)(Sl 2 (Q)(z)-Sr(Q)(z))dv(z) 
< — k Y. h{\T{2-)\, |Q 12 |, |Q 21 |, . . . , \Q ab - Q a % . . . , |Q 12 |) 

3=0 

where the difference occurs in the j'th spot and ab = 12 if j is even or 21 if j is odd. 
Thus, Lemma El implies that if \\QW2 < M and \\QW2 < M, then for each e > 0, there 
is a constant C f so that 



\S 12 (Q) - S 12 (Q)\\ 2 <C e \\Q- Q\\ 2 M 2 £(2fc + 1)(1 + e) 



2A- 



1 



fc=0 



M 2(fc-1)_ 



Since we have M < a/2, we have the Lipschitz continuity also. 

The inverse map is handled similarly so we only give a sketch of the argument. 
The map which takes Q to S is given by 

Q( x ) = -J f Tm{x,z)dfi{z) (18) 

7T JC 

where the map T depends on S and was defined in (JHJ). We have that m may be 
represented as 

m= (I -CT)-\1) 



where C is the Cauchy transform 

Cf{z) 



fW 



7T JC Z — W 



dfi(w) 
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acting on matrix valued functions. As before, we can write m as the series, 



m=£(CT)'(l) (19) 

3=0 

Substituting (fTH|) into (fTH|) gives a series representation for Q, Q(x) = - J2kLoQk(x) 
where the kth term is given by 

J_ r s 12 (z )s 21 ( Zl )...s 12 (z 2k ) 

7r2fc J C 2k+i (z Q - z{){z\ - z 2 ) ... (hk-l ~ z 2k) 

xa 1 (x,z -z 1 + ...- z 2k -i + z 2k ) d/j,(z , . . . z 2k ). 

Arguing as before, we can show that for each e > 0, there is constant C e so that 

\\Q2 k h < cS^-^\\s\\l k+1 

and also obtain that Q depends continuously on S. 

The statements that Q o S = Id and S o Q = Id follow from continuity and the 
corresponding results for nice potentials as proved by Sung [TBI EH HE] ■ However, 
one must replace Lemma 2.3 of jTHj. The uniqueness of m needed to carry Sung's 
arguments holds for ||Q||2 < \/2, thanks to the estimate for at (3/4) in LemmaQ(see 
also l)17p). Also, the extension of the non-linear Plancherel identity to functions in 
L 2 is immediate from the result for smooth functions and the continuity of the map 
Q - S. i 

We give an application of our estimate. This application depends on the fact that 
we have shown that the scattering map and its inverse map a neighborhood of in 
L 2 into L 2 . Since the domain and range lie in the same space, L 2 , it is particularly 
simple to construct the solution. 

Recall that the Davey-Stewartson II equation as studied by Beals and Coifman is 
the evolution equation in two space dimensions given by 

q t = iq x i x 2 - Airq 

Tx 1 x 1 Tx 2 x 2 (\q\ )x 1 x 2 

If we set 

q(x,t) \ 
q(x,t) ) 

and apply the scattering map to obtain S((-,t) = S(Q)(-,t), then the function S(t) 
satisfies 

S t = (4zz 1 z 2 )5. 

Thus the solution of the Davey-Stewartson II equation is the upper-right entry of the 
matrix 

Q(-,t) = Q(exp(4^^ 2 )5(Q(-, 0))). (20) 



Q 12 (x) = 



Q(x,t) = 
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This is established for nice potentials Q by Beals and Coifman [3], see also Sung 
for a detailed argument El HM ■ O ur contribution is the observation that the 
maps Q and S are continuous in L 2 which immediately implies the following result. 
See Ghidaghlia and Saut jH] and Linares and Ponce ^Hj for additional results on the 
Davey-Stewartson systems. 

Corollary 5 The scattering solution of the Davey-Stewartson II system as defined 
in gives a solution which depends continuously on the initial data in L 2 . If the 
initial data qj(-,0) satisfies \\qj(-,0)\\2 < 1 an d say qj(x,0) G S(C) for j = 1 and 2, 
then the solutions qj (x, t) satisfy 

\\qi(;t)-q 2 (;t)\\2<C\\ qi (;0)-q 2 (;0)\\2- 
We close with three questions related to the above results. 

1. Can we remove the smallness restriction from our Theorem It is known that 
if one of Q = ±Q* holds and Q G L 2 , then the operator / — G Z Q is invertible. 
This fact does not seem to be in the literature. To see the invertibility of I—G Z Q 
on we observe that arguments similar to those of Nachman fH\ Lemma 
4.2] show that the map / — > G z Qf is compact on L 4//3 . The injectivity of the 
map /—>•/ — G z Qf when Q = ±Q* follows from Corollary 3.11 of This 
Corollary is an observation of Nachman's. Hence, the Fredholm theory gives 
the claim. 

2. Can we construct the Jost solutions m when Q is only assumed to be in L 2 ? 

3. The expression (|2Uj) is defined, if Q is just in L 2 and small. In what sense does 
this expression solve the Davey-Stewartson equations? 
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